Low Energy Neutron Production by Inverse f3 decay in Metallic Hydride Surfaces 
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It has been recently argued that inverse-/? nuclear transmutations might occur at an impressively 
high rate in a thin layer at the metallic hydride surface under specific conditions. In this note we 
present a calculation of the transmutation rate which shows that there is little room for such a 
remarkable effect. 

PACS: 24.60.-k, 23.20.Nx 



Introduction and main results. In a number of papers [lH3j, A. Widom and collaborators consider the 
intriguing possibility that low energy neutrons are produced on metallic hydride surfaces by the capture of 
electrons by protons and the subsequent inverse-/3 reaction: 

ep — > nv e (1) 

This process cannot, of course, take place in a Hydrogen atom in vacuum, it would have a negative Q-value 
(Q = Mp+m e — M n < 0). The authors advance the hypothesis that mass renormalization from the interaction of 
electrons in atoms with some electromagnetic external field (e.g. due to a laser source) may increase the effective 
electron mass, so as to make it possible the reaction (Q]). The tilde sign on the electron symbol underscores that 
the electron bound to the proton is a 'dressed' one, with an effective mass denoted by: 

m* = f3m e (2) 

A central element of the analysis of this very intriguing effect is the decay rate of the bound ep state which 
is given in Ref. [|[ as 0: 

T WL (ep -+ n Ve ) = {GF o m } )2 (l + 3A 2 )m e (/3 - /3 ) 2 (3) 

Natural units are used throughout, we approximate to unity the cosine of the Cabibbo angle, A = gA/gv ~ 1-25 
is the axial to vector coupling in neutron /3 decay. 
The threshold value for process ([TJ) to occur at all is: 

Po = — - M " « 2.59 (4) 

m e 

With the above numerical values and j3 sa 20, the authors find: 

r wL -> nv e ) « (1.86 x 1(T 3 ) X (fi - A,) 2 Hz « 0.56 Hz (5) 

What is surprising in formula ((SJ is the absence of the factor |-!/>(0)| 2 , which gives the probability for finding 
the electron and the proton at the same point Q. This factor is a mandatory consequence of the local Fermi 
lagrangian density: 



£(x) = ~j= [ n (x)ln( l ~ ^5)p(x)] [^OOTmC 1 ~ 75)e(x)] (6) 

With the electron and proton fields computed at the same point, the factor |^(0)| 2 indicates that it must be 
difficult for the lagrangian to induce transitions from an initial state with, say, the proton localized in our lab 
and the electron on Mars. We present a calculation of the rate of ([T]) done in two independent ways. 

One is to consider the decay of the ep bound states with spin and 1, relating the amplitude to the basic 
lagrangian © in a way analogous to the quark model calculation of pion decay made time ago by Van Royen 
and Weisskopf (VRW) 0]. 



The second is to use the time-honored formula: 



r = \r/>(o)\ 2 va (7) 

where a is the unpolarized cross section for process ([1]). This is the basic formula in the theory of electron 
K-capture, also used [B[ to derive the decay rate of positronium from the annihilation cross section of an 
electron-positron pair into photons. 
Either ways we find, reassuringly: 

T(ep -> nv e ) = |^(0) | 2 x i- (G F m e f (l + 3A 2 ) (J3 - f3 ) 2 (8) 

Defining: 

l^(0)| 2 = 1/R 3 (9) 

the result in ([3]) is consistent with Eq. (|SJ) for: 

R l/m e ~ 0.4 x 10~ 2 A (10) 

that is the electron should be confined within its Compton radius, which is completely unrealistic. 
We further set 

7r(a*) c> 7ra J 

with (3 defined in ([2]), a the Bohr radius, a = (arn e ) _1 ~ 0.54 A and a the fine structure constant, and find, 
finally: 



T{ep -> ni/ e ) = ^2 ™') 2 ™ e (1 + 3A 2 ) /? 3 (/3 - f3 ) 2 (12) 

Numerically: 

r(Sp-> ra/ e @/3 = 20) = 1.8 x 10 -3 Hz (13) 

while, for a more moderate value, j3 « 2/3o: 

r(ep-> nv e @P = 5.2) = 6.9 x 10~ 7 Hz (14) 

The numerical value proposed in PHI], Eq. (J5J, is obtained for /3 = 61 (i.e. m* » 30 MeV). 

In the following, we give the details of the calculation of T. In the end, we comment on the high values of j3 
considered in [1-y]. 

Decay rates of S-wave bound states. Neutron emission can be described as arising from the decay of the 
two, S-wave ground states, Hq and H\, with total spin equal to zero and one, respectively. We may describe 
this decay with two phenomenological parameters, /o,i according to: 

Cphen = -7= [f n(x)(l +"f 5 )v c H a + /in(x)7i(l + ~/ 5 )v c H[] (15) 

The Fermi constant has been inserted for convenience, Hq and H\ are field operators describing the annihilation 
of either H state and we describe the creation of the neutrino with the antineutrino field, v c . The situation is 
entirely similar to the decay 7r~ — > ii~v c , described by the phenomenological parameter / w [To| . 

To connect to the VRW formulae, it is convenient to rewrite the Fermi lagrangian so as to have the annihilation 
operators for the electron and the proton in the same Dirac bilinear (similarly to the quark-antiquark fields in 
pion decay). This is done by first introducing the positron and antineutrino fields in © and then making a 
Fierz transformation. To wit [Tlj : 

[^Ob^ 1 - 75)e(x)] [n(x)7 M (l - Xj 5 )p{x)} = - [e c (x)^(l + ^ 5 )v c (x)] [fi(x)^(l - X-f 5 )p(x)} = (16) 

! " A [e c (x)(l - j 5 )p(x)j [n(x)(l + j 5 )v c (x)} - (^-) [e c {x)f (1 + l5 )p(x)] [n(x)yi(l + 1s)v c (x)\ 
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Following VRW we define the correctly normalized H state according to: 

\H) = (2^) 3 / 2 / d 3 p /(p)0(r, S )K)t(p)K)t(-p)|O) (17) 



r and s are spin indices, </>(r, s) are the Clebsch-Gordan coefficients appropriate to the spin of H, (a P:6 )^ are the 
proton and electron creation operators and f(p) the momentum space wave function, with: 

Jd 3 P \f(p)\ 2 = l (18) 

The parameter f is obtained by comparing (fT5|) with (fl~7|). e.g.: 

fo(0\H \H ) = (0| [e c (0)(l - 7 5>(0)] \H ) (19) 



Expanding the fields in normal modes, spinors and gamma matrices combine with the appropriate Clebsches, 
and we obtain (the y/2 factor comes from the Clebsch, the overall signs of fo.i are irrelevant): 



with the x-space wave function given by: 



= (^72 / *P «-**7(p) (21) 



/i = -(^) V2^(0) (22) 



Similarly: 

Direct calculation of the spin-averaged rates gives: 

r ( H„^ „ + „) . 4 g££ <li^l! W o)|^ - A )> 

r(ff, -, n + ,) = 4 <2^£>! <l^l!| f (0)l' tf - A )» (23) 
and the inclusive rate is obtained from: 

T(gp -> ni/ e ) - i [T(flo -> n + i/) + 3r(H x -> n + i/)] = |^(0)| 2 £*£fhL(i + 3A 2 )(/3 - /3 ) 2 (24) 

4 Z7T 

as anticipated in j8]). 

Inclusive rate from the unpolarized cross section. The cross section for e~ + p ^ ^ + n is @: 

(G f m e ) 2 , ~ (/3-/3 ) 2 
cr(e +p^^ + n) = (1 + 3A) (25) 

The transmutation rate thus obtained via Eq. ([7]) confirms the result in (|8|). 

TTie value of (3. Values of /3 of the order or even larger than twenty are certainly unusual in condensed 
matter physics, expecially for bound electrons. An estimate of (3 is given in [H-Q 



where A is given in terms of the plasma frequency Q p of the protons 

M 



A=^ (MQl)a 2 (27) 



3 



and y \u\ is the r.m.s. displacement of the protons. 




3]: 



Current values of fl v are estimated to be of order fl 



■v 



0.1 eV 0. With this value of fl p and J \u\ /a 2 




4.2 0- 



■p 



0.1) = 1.01 



(28) 



The value Q p = 0.8 eV used in flji3| leads to: 



■v 



0.8) = 1.8 



(29) 



considerably lower than j3 ~ 20 and in any case below threshold for nuclear transmutation to occur. 

Conclusions. A correct calculation gives a neutron production rate from (pj about 300 times smaller than 
what estimated in [lJ-Q , for the value of the mass renormalization factor j3 ~ 20 considered there. In turn, 
it is questionable that values of j3 can be realized, in particular for bound electrons, so large as to give rise to 
useful nuclear transmutation rates. A more detailed analysis of the attainable values of f3 is needed to obtain 
more definite conclusions on this interesting phenomenon, should it exist at all. 
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